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Recent developments in the study of primordial black holes (PBHs) will be reviewed, with partic- 
ular emphasis on their formation and evaporation. PBHs could provide a unique probe of the early 
Universe, gravitational collapse, high energy physics and quantum gravity. Indeed their study may 
place interesting constraints on the physics relevant to these areas even if they never formed. In 
the "early Universe" context, particularly useful constraints can be placed on inflationary scenarios, 
especially if evaporating PBHs leave stable Planck-mass relicts. In the "gravitational collapse" con- 
text, the existence of PBHs could provide a unique test of the sort of critical phenomena discovered 
in recent numerical calculations. In the "high energy physics" context, information may come from 
gamma-ray bursts (if a subset of these are generated by PBH explosions) or from cosmic rays (if 
some of these derive from evaporating PBHs). In the "quantum gravity" context, the formation 
and evaporation of small black holes could lead to observable signatures in cosmic ray events and 
accelerator experiments, providing there are extra dimensions and providing the quantum gravity 
scale is around a TeV. 



I. INTRODUCTION 

Black holes with a wide range of masses could have 
formed in the early Universe as a result of the grea t 
compression associated with the Big Bang [7(3, Il43| . 
A comparison of the cosmological density at a time t 
after the Big Bang with the density associated with 
a black hole of mass M shows that PBHs would have 
of order the particle horizon mass at their formation 
epoch: 



M H {t) 



G 



10 



15 



t 



10 



-23 



9- 



(1) 



PBHs could thus span an enormous mass range: those 
formed at the Planck time (10 _43 s) would have the 
Planck mass (10 _5 g), whereas those formed at 1 s 
would be as large as 10 5 M©, comparable to the mass 
of the holes thought to reside in galactic nuclei. By 
contrast, black holes forming at the present epoch 
could never be smaller than about 1M©. 

The realization that PBHs might be small prompted 
Hawking to study their quantum properties. This led 
to his famous discovery [7l| that black holes radiate 
thermally with a temperature 
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so they evaporate on a timescale 
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Only black holes smaller than 10 15 g would have evapo- 
rated by the present epoch, so eqn (1) implies that this 
effect could be important only for black holes which 
formed before 10 _23 s. 

Hawking's result was a tremendous conceptual ad- 
vance, since it linked three previously disparate areas 



of physics - quantum theory, general relativity and 
thermodynamics. However, at first sight it was bad 
news for PBH enthusiasts. For since PBHs with a 
mass of 10 15 g would be producing photons with en- 
ergy of order 100 MeV at the present epoch, the ob- 
servational limit on the 7-ray background intensity at 
100 MeV immediately implied that their de nsity could 
not exceed 10 -8 times the critical density |l22j |. Not 
only did this render PBHs unlikely dark matter can- 
didates, it also implied that there was little chance of 
dete cting black hole explosions at the present epoch 

EH. 



Nevertheless, it was soon realized that the 7-ray 
background results did not preclude PBHs playing 
other important cosmological roles [27) and some of 
these will be discussed in this review. Indeed the dis- 
covery of a PBH could provide a unique probe of at 
least four areas of physics: the early Universe; grav- 
itational collapse; high energy physics; and quantum 
gravity. The first topic is relevant because studying 
PBH formation and evaporation can impose impor- 
tant constraints on primordial inhomogeneities and 
cosmological phase transitions. The second topic re- 
lates to recent developments in the study of "critical 
phenomena" and the issue of whether PBHs are vi- 
able dark matter candidates. The third topic arises 
because PBH evaporations could contribute to cos- 
mic rays, whose energy distribution would then give 
significant information about the high energy physics 
involved in the final explosive phase of black hole 
evaporation. The fourth topic arises because it has 
been suggested that quantum gravity effects could ap- 
pear at the TeV scale and this leads to the intrigu- 
ing possibility that small black holes could be gener- 
ated in accelerators experiments or cosmic ray events, 
with striking observational consequences. Although 
such black holes are not technically "primordial" , this 
possibility would have radical implications for PBHs 
themselves. 
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II. PBHS AS A PROBE OF PRIMORDIAL 
INHOMOGENEITIES 

One of the most important reasons for studying 
PBHs is that it enables one to place limits on the 
spectrum of density fluctuations in the early Universe. 
This is because, if the PBHs form directly from den- 
sity perturbations, the fraction of regions undergo- 
ing collapse at any epoch is determined by the root- 
mean-square amplitude e of the fluctuations entering 
the horizon at that epoch and the equation of state 
p = 7/5 (0 < 7 < 1). One usually expects a radiation 
equation of state (7 = 1/3) in the early Universe but 
it may have deviated from this in some periods. 



A. Simplistic Analysis 

Early calculations assumed that the overdense re- 
gion which evolves to a PBH is spherically symmetric 
and part of a closed Friedmann model. In order to 
collapse against the pressure, such a region must be 
larger than the Jeans length at maximum expansion 
and this is just ^7 times the horizon size. On the 
other hand, it cannot be larger than the horizon size, 
else it would form a separate closed universe and not 
be part of our Universe p9j. 

This has two important implications. Firstly, PBHs 
forming at time t after the Big Bang should have of 
order the horizon mass given by eqn (1). Secondly, for 
a region destined to collapse to a PBH, one requires 
the fractional overdensity at the horizon epoch 5 to 
exceed 7. Providing the density fluctuations have a 
Gaussian distribution and are spherically symmetric, 
one can infer that the fraction of regions of mass M 
which collapse is [2^ 



P(M) ~ e(M) exp 



2e{Mf 



(4) 



where e(M) is the value of e when the horizon mass is 
M. The PBHs can have an extended mass spectrum 
only if the fluctuations are scale-invariant (i.e. with 
e independent of M). In this case, the PBH mass 
distribution is given by [26j 

dn/dM = (a- 2)(M/M,)- a M- 2 n PBH p clit (5) 

where M* as 10 15 g is the current lower cut-off in the 
mass spectrum due to evaporations, Qpbh is the total 
density of the PBHs in units of the critical density 
(which itself depends on /?) and the exponent a is 
determined by the equation of state: 
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a = 5/2 if one has a radiation equation of state. This 
means that the density of PBHs larger than M falls 



off as M 1 / 2 , so most of the PBH density is contained 
in the smallest ones. 

Many scenarios for the cosmological density fluctu- 
ations predict that e is at least approximately scale- 
invariant but the sensitive dependence of /3 on e means 
that even tiny deviations from scale-invariance can be 
important. If e(M) decreases with increasing M, then 
the spectrum falls off exponentially and most of the 
PBH density is contained in the smallest ones. If e(M) 
increases with increasing M , the spectrum rises expo- 
nentially and - if PBHs were to form at all - they could 
only do so at large scales. However, the microwave 
background anisotropics would then be larger than 
observed, so this possibility can be rejected. 

The current density parameter ^pbh associated 
with PBHs which form at a redshift z or time t is 
related to by 
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where flu as 10 4 is the density parameter of the mi- 
crowave background and we have used eqn (1). The 
(1 + z) factor arises because the radiation density 
scales as (1 + z) 4 , whereas the PBH density scales 
as (1 + z) 3 . Any limit on fipBH therefore places a 
constraint on 0{M) and the constraints are summa- 
rized in Fig. ^ which is taken from Carr et al. |32^ |. 
The constraint for non-evaporating mass ranges above 
10 15 g comes from requiring SIpbh < 1 but stronger 
constraints are associated with PBHs smaller than 
this since they would have evaporated by now. The 
strongest one is the 7-ray limit associated wit h th e 
10 15 g PBHs evaporating at the present epoch .122]. 
Other ones are associated with the generation of en- 
tropy and modification s to the cosmological produc- 
tion of light elements |120| . The constraints below 
10 6 g are based on the (uncertain) assumption that 
evaporating PBHs leave stable Planck mass relics, an 
issue which is discussed in detail in Section 4 . Other 
constraints, not shown here, are associated with grav- 
itin o pro duction [9l| . reionization |75j| . thermodynam- 
ics |l0lj ] and the holographic principle . 

The constraints on 0(M) can be converted into con- 
straints on e(M) using eqn (4) and these are shown 
in Fig. |21 Also shown here are the (non-PBH) con- 
straints associated with the spectral distortions in the 
cosmic microwave background induced by the dissipa- 
tion of intermediate scale density perturbations and 
the COBE quadrupole measurement. This shows that 
one needs the fluctuation amplitude to decrease with 
increasing scale in order to produce PBHs and the 
lines corresponding to various slopes in the e(M) re- 
lationship are also shown in Fig. [3 
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FIG. 1: Constraints on /3(M) 
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FIG. 2: Constraints on e(M) 



B. Refinements of Simplistic Analysis 

The criterion for PBH formation given above is 
rather simplistic and needs to be tested with de- 
tailed numerical calculations. The first hydrodynam- 
ical studies of P BH formation were carried out by 
Nadezhin et al. |l!5l |. These roughly confirmed the 
criterion S > 7 for PBH formation, although the PBHs 
could be be somewhat smaller than the horizon. In 
recent years several groups have carried out more de- 
tailed hydrodynamical calculations and these have re- 
fined the S > 7 criterion and hence the estimat e for 
(3{M) given by eqn (4). Niemeyer & Jedamzik 
find that one needs S > 0.7 rather than 5 > 0.3 to 
ensure PBH formation and they also find that there is 
little accretion after PBH formatio n, as expected the- 
oretically pflj . Shibata & Sasaki |l32j reach similar 
conclusions. 

It should be stressed that the description of fluc- 
tuations beyond the horizon is somewhat problematic 
and it is clearer to use a gauge-invariant description 
whic h inv olves the total energy or metric perturba- 
tion |l32fl . Also the derivation of the mass spectrum 



given by eqn (4) is based on Press-Schechter theory 
and it is more satisfactory to use peaks theory. Both 
these points have been considered by Green et al. |fi3 |. 
They find that that the critical value for the density 
contrast is around 0.3, which (ironically) is close to 
the value advocated 30 years ago! 

Another refinement of the simplistic analysis which 
underlies eqn (4) concerns the assumption that the 
fluctuations have a Gaussian distribution. Bullock & 
Primack [24| and Ivanov [8l| have pointed out that 
this may not apply if the fluctuations derive from an 
inflationary period (discussed in more detail later). 
So long as the fluctuations are small (8<p/4> <C 1), as 
certainly applies on a galactic scale, this assumption 
is valid. However, for PBH formation one requires 
8<j)/(f) ~ 1, and, in this case, the coupling of different 
Fourier modes destroys the Gaussianity. Their analy- 
sis suggests that (3{M) can be very different from the 
value indicated by eqn (4) but it still depends very 
sensitively on e. 



C. PBHs and Critical Collapse 

A particularly interesting development has been the 
application of "critical phenomena" to PBH forma- 
tion. Studies of the collapse of various types of spher- 
ically symmetric matter fields have shown that there 
is always a critical solution which separates those con- 
figurations which form a black hole from those which 
disperse to an asymptotically flat state. The config- 
urations are described by some index p and, as the 
critical index p c is approached, the black hole mass 
is found to scale as (p — p c ) n for some exponent 77. 
This effect was first discovered for scalar fields 
but subsequently demonstrated for radiation |49j ] and 
the n mo re general fluids with equation of state p = 

mini. 

In all these studies the spacetime was assumed to be 
asym ptotically flat. However, Niemeyer & Jedamzik 
|I17| have applied the same idea to study black hole 
formation in asymptotically Friedmann models and 
have found similar results. For a variety of initial den- 
sity perturbation profiles, they find that the relation- 
ship between the PBH mass and the the horizon-scale 
density perturbation has the form 



M = KM H {5 — S c ) v 



(8) 



where Mu is the horizon mass and the constants are 
in the range 0.34 < -q < 0.37, 2.4 < K < 11.9 
and 0.67 < S c < 0.71 for the various configura- 
tions. Since M — ► as S — > 6 C , this suggests that 
PBHs may be much smaller than the particle hori- 
zon at formation and it also modifies the mass spec- 
trum [H IHJ. However, recently Miller et al. 
|114| have found that the critical overdensity lies in 
the lower range 0.43 < S c < 0.47 if one only allows 
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growing modes at decoupling (which is more plausible 
if the fluctuations derive from inflation). They also 
find that the exponent n is modified if there is a cos- 
mological constant. Hawke & Stewart [t^j claim that 
the formation of shocks prevents black holes forming 
on scales belo w 10~ 4 of the horizon mass but this has 
been disputed |l!4j . 



III. PBHS AS A PROBE OF THE EARLY 
UNIVERSE 

Many phase transitions could occur in the early 
Universe which lead to PBH formation. In some of 
these one require pre-existing density fluctuations but 
in others the PBHs form spontaneously even if the 
Universe starts off perfectly smooth. In the latter 
case, f3(M) depends not on e(M) but on some other 
cosmological parameter. 



C. Bubble Collisions 

Bubbles of broken symmetry might arise at any 
spontaneously broken symmetry epoch and various 
people have suggested that PBHs could form as a re- 
sult of bubble collisions |43l u4 . l99j . However, this 
happens only if the bubble formation rate per Hubble 
volume is finely tuned: if it is much larger than the 
Hubble rate, the entire Universe undergoes the phase 
transition immediately and there is not time to form 
black holes; if it is much less than the Hubble rate, 
the bubbles are very rare and never collide. The holes 
should have a mass of order the horizon mass at the 
phase transition, so PBHs forming at the GUT epoch 
would have a mass of 10 3 g, those forming at the elec- 
troweak unification epoch would have a mass of 10 28 g, 
and those forming at the QCD (quark-hadron) phase 
transition would have mass of around 1M@. 

D. PBHs and Inflation 



A. Soft Equation of State 

Some phase transitions can lead to the equation of 
state becoming soft (7 -C 1) for a while. For example, 
the pressure may be reduced if the Universe's mass 
is ever channelled into particles which are massive 
enough to be non-relativistic. In such cases, the effect 
of pressure in stopping collapse is unimportant and 
the probability of PBH formation just depends upon 
the fraction of regions which are sufficiently spherical 
to undergo collapse |90| . For a given spectrum of pri- 
mordial fluctuations, this means that there may just 
be a narrow mass range - associated with the period 
of the soft equation of state - in which the PBHs form. 



B. Collapse of Cosmic Loops 

In the cosmic string scenario, one expects some 
strings to self-intersect and form cosmic loops. A typ- 
ical loop will be larger than its Schwarzschild radius 
by the factor (G^i) -1 , where fi is the string mass per 
unit length. If strings play a role in generating large- 
scale structure, G\i must be of order 10~ 6 . However , 
as discussed by many authors [H IH El Efll Il2^ |. 
there is always a small probability that a cosmic loop 
will get into a configuration in which every dimension 
lies within its Schwarzschild radius. This probability 
determines the collapse fraction p and depends upon 
both /1 and the string correlation scale. Note that eqn 
(5) still applies since the holes are forming with equal 
probability at every epoch. 



Inflation has two important consequences for PBHs. 
On the one hand, any PBHs formed before the end of 
inflation will be diluted to a negligible density. In- 
flation thus imposes a lower limit on the PBH mass 
spectrum: 

M > M min = M P {T RH /T Pl )- 2 (9) 

where Trh is the reheat temperature and Tp « 
10 19 GeV is the Planck temperature. The CMB 
quadrupole measurement implies Trh ~ 10 16 GeV, so 
-Mmin certainly exceeds 1 g. On the other hand, infla- 
tion will itself generate fluctuations and these may suf- 
fice to produce PBHs after reheating. If the inflaton 
potential is V(<fi), then the horizon-scale fluctuations 
for a mass-scale M are 

/ y3/2 \ 

Cma {MW) H <10) 

where a prime denotes d/d<j> and the right-hand side 
is evaluated for the value of <j> when the mass-scale 
M falls within the horizon. In the standard chaotic 
inflationary scenario, one makes the "slow-roll" and 
"friction-dominated" asumptions: 

£ = {MpV'/Vf < 1, n = MpV"/V < 1. (11) 

Usually the exponent n characterizing the power spec- 
trum of the fluctuations, |<5fc| 2 « k n , is very close to 
but slightly below 1: 

n=l+4£- 277 wl. (12) 

bmcc c scales as M^ 1 "™^ 4 , this means that the fluc- 
tuations are slightly increasing with scale. The nor- 
malization required to explain galaxy formation (e w 
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10 -5 ) would then preclude the formation of PBHs 
on a smaller scale. If PBH formation is to occur, 
one needs the fluctuations to decrease with increasing 
mass (n > 1) and, from eqn (13), this is only possible 
if the scalar field is accelerating sufficiently fast that 

V"/V > (1/2){V'/V) 2 . (13) 

This condition is certainly satisfied in some scenarios 
[30| and, if it is, eqn (4) implies that the PBH density 
will be dominated by the ones forming immediately 
after reheating. 

Since each value of n corresponds to a straight line 
in Fig|Sl any particular value for the reheat time ti 
corresponds to an upper limit on n. This limit is indi- 
cated in FigOU which is taken from [23, apart from a 
correction pointed out by Green & Liddle |H3|- Sim- 
ilar constraints have been obtained by several other 
people 0, . The figure also shows how the con- 
straint on n is strengthened if the reheating at the end 
of inflation is sufficiently slow for there to be a dust- 
like phase [62 . However, it should be stressed that 
not all inflationary scenarios predict that the spectral 
index should be constant. Hodges & Blumcnthal |78| 
have pointed out that one can get any spectrum for 
the fluctuations whatsoever by suitably choosing the 
form of V{4>). For example, eqn (10) suggests that 
one can get a spike in the spectrum by flattening the 
potential over some mass range (since the fluctuation 
diverges when V goes to 0). This idea was exploited 
by Ivanov et al. [82|. who fine-tuned the position of 
the spike so that it corresponds to the mass-scale asso- 
ciated with microlcnsing events observed in the Large 
Magellanic Cloud. 




FIG. 3: Constraints on spectral index n in terms of reheat 
time ti 

It should be noted that the relationship between the 
variance of the mass fluctuations relevant for PBHs 
and the present day horizon-scale density fluctuations 
in the inflation scenario is not trivial. For scale-free 
fluctuations, there is a simple relationship between 
them but it is more complicated otherwise. Blais et 



al. [l9j find that the mass fluctuations are reduced by 
34% for a spec tral index in the range 1 < n < 1.3, 
while Polarski |l23j | finds a further 15% reduction if 
there is a cosmological constant. 

Even if PBHs never actually formed as a result of 
inflation, studying them places important constraints 
on the many types of inflationary scenarios. Besides 
the chaotic scenario discussed above, there are al so th e 
variants of inflation described as sup ernatural [126], 
supersym metric [57|. hybrid l87l. m ultiple |140| . 
oscill ating |134| , preheating |l3l l47l l52l lo3l ] , running 
mass |l00j and saddle • There are also scen arios in 
which the inflaton serves as dark matter |l02j . PBH 
formation has been studied in all of these models. 



E. Varying Gravitational Constant 

The PBH constraints would be severely modified if 
the value of the gravitational "constant" G was dif- 
ferent at early times. The simplest varying-G model 
is Brans-Dicke (BD) theory J2.2] , in which G is associ- 
ated with a scalar field <j) and the deviations from gen- 
eral relativity are specified by a parameter cu. A vari- 
ety of astrophysical tests currently require |u>| > 500, 
which implies that the deviations can only ever be 
small |l37j . Howev er, there e xist generalized scalar- 
tensor theories ll 19i Il36| in which uj is itself a 
function of <fi and these lead to a considerably broader 
range of variations in G. In particular, it permits u> 
to be small at early times (allowing noticeable varia- 
tions of G then) even if it is large today. In the last 
decade interest in such theories has been revitalized as 
a result of early Universe studies. Extended inflation 
explicitly requires a model in which G varies [9!| and, 
in higher dimensional Kaluza-Klein-type cosmologies, 
the variation in the sizes o f the extra dimensions also 
naturally leads to this [53l lllC| . 

The consequences of the cosmological variation of 
G for PBH evaporation depend upon how the value 
of G near the black hole evolves. Barrow ^(| intro- 
duces two possibilities: in scenario A, G everywhere 
maintains the background cosmological value (so <p is 
homogeneous); in scenario B, it preserves the value 
it had at the formation epoch near the black hole 
even though it evolves at large distances (so 4> be- 
comes inhomogeneous). On the assumption that a 
PBH of mass M has a temperature given by eqn (2) 
with G = G(t) in scenario A and G = G(M) in sce- 
nario B, Barrow & Carr [ljj calculate how the evap- 
oration constraints summarized in Fig.l are modified 
for a wide range of varying-G models. The question 
of whether scenario A or scenario B is more plausible 
has been considered in several papers |28t |83j but re- 
cent studies of the collapse of a scalar field to a PBH 
suggest that gravitational memory is unlikely [6?], [6^ . 
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IV. PBHS AND DARK MATTER 



B. PBHs and Supermassive Black Holes 



Roughly 30% of the total density of the Universe is 
now thought to be in the form of "cold dark mat- 
ter". Recently there has been a lot of interest in 
whether PBHs could provide this, since those larger 
than 10 15 g would not have evaporated yet and would 
certainly be massive enough to be dynamically "cold" . 
They might alsoplay a rol e in the formation of large- 
scale structure 0, |23> 1.109]. However, there are some 
ranges in which this is excluded. For example, fem- 
tolensing of gamma-ray bursts by PBHs precludes 
those in the mass range 10 17 — 10 20 g from having 
a critical density and searches for microlensing of 
stars in the Large Magellanic Cloud, while allowing 
a tenth-critical-density at around a solar mass, ex- 
clude 10 26 - 10 34 g PBHs 0. However, there are no 
constraints in the intermediate (sublunar) mass range 

10 20 - 10 26 g m. 



A. PBH Formation at Phase Transitions 

One possibility is that PBHs with a mass of around 
IMq could have formed at the quark-hadron phase 
transition at 10 _5 s because of a temporary softening 
of the equation of state then • Such PBHs would 
naturally have the sort of mass required to explain the 
MACHO microlensing results If the QCD phase 
transition is assumed to be 1st order, then hydrody- 
namical calculations show that the value of 6 required 
for PBH formation is indeed reduced below the value 
which pertains in the radiation case 85]. This means 
that PBH formation will be strongly enhanced at the 
QCD epoch, with the mass distribution peaking at 
around the horizon mass then. Another possibility is 
that PBHs with a mass of around 10 _7 M Q could form 
in TeV quantum gravity scenarios |79j . 

One of the interesting implications of these scenar- 
ios is the possible exis tence of a halo population of 
binary black holes |ll6j |. With a full halo of such ob- 
jects, there could be a huge number of binaries in- 
side 50 kpc and some of these could be coalescing due 
to g ravitational radiation losses at the present epoch 
20]. If the associated gravitational waves were de- 
tected, it would provide a unique probe of the halo 
distribution [8(j. Gravity waves from binary PBHs 
would be detectable down to 10 _5 Mq using VIRGO, 
10~ 7 M Q using EURO and 10~ n M Q using LISA [fjj. 
Note that LISA could also detect isolated PBHs by 
measuring the gravitati onal impulse induced by any 
nearby passing ones P, Il3l| . However, this method 
would not work below 10 14 g (because the effect would 
be hidden by the Moon) or above 10 20 g (because the 
encounters would be too rare). 



Several people have suggested that the 10 6 — 1O 8 M 
supermassive black holes thought to resid e in g alactic 
nuclei could be of primordial origin 46. 86l lT27f . How- 
ever, since no PBHs are likely to form after 1 s, cor- 
responding to a maximum formation mass of 10 5 Mq, 
this requires a large amount of accretion. Bean and 
Magueijo have proposed [14( that PBHs may accrete 
from the quintessence field which is invoked to explain 
the acceleration of the Universe. However, th ey in voke 
a Newtonian formula for the accretion rate jl43j ] and 
this is known to be questionable [29j. Recent studies 
of the accretion of a scalar field by a PBH also indicate 
that this is unlikely |ol|. 



C. Planck Mass Relics 

Some people have speculated that black hole evap- 
oration could cease once the hole gets close to the 
Planck mass 0, ^J. For example, in the stan- 
dard Kaluza-Klein picture, extra dimensions are as- 
sumed to be compactified on the scale of the Planck 
length. This means that the influence of these extra 
dimensions becomes important at the energy scale of 
10 19 GeV on which quantum gravity effects become 
significant. Such effects could influence black hole 
evaporation and could conceivably result in evapo- 
ration ceasing at the Planck mass. Various non- 
quantum-gravitational effects (such as higher order 
corrections to the gravitational Lagrangian or string 
effects) could also lead to stable relics [33 but the relic 
mass is usually close to the Planck scale. 

Another possibility, as argued by Chen & Adler [34[ , 
is that stable relics could arise if one invokes a "gener- 
alized uncertainty principle" . This replaces the usual 
uncertainty principle with one of the form 



Ax > -j^ + l P ^-, 



n 

Ap 



(14) 



where the second term is supposed to account for self- 
gravity effects. This means that the black hole tem- 
perature becomes 



Tbh = 



Mc 2 
ink 



1 - 



Ml 
M 2 



(15) 



This reduces to the standard Hawking form for M S> 
Mp but it remains finite instead of diverging at the 
Planck mass itself. 

Whatever the cause of their stability, Planck mass 
relics would provide a possible cold dark matter candi- 
date |l03j ] . Indeed this leads to the "relics" constraints 
indicated in Fig.l, Fig. 2 and Fig. 3. In particular, such 
relics could be left over from inflation 

[HOI. If the 
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relics have a mass ttMp, then the requirement that 
they have less than the critical density implies 

(3{M) < IQ^n^iM/Mpf' 2 (16) 

for the mass range 

(Trh/Tp)- 2 M pi <M< 10 n K 2/5 A/ P . (17) 

The upper mass limit arises because PBHs larger than 
this dominate the total density before they evapo- 
rate. Producing a critical density of relics obviously 
requires fine-tuning of the index n. Also one needs 
n sal. 3, which is barely compatible with the WMAP 
results 0, ■ Nevertheless, this is possible in princi- 
ple. In particular, Chen |33j has argued that hybrid 
inflation could produce relics from 10 5 g PBHs formed 
at l(T 32 s. 



V. PBHS AS A PROBE OF HIGH ENERGY 
PHYSICS 

A black hole of mass M will emit particles like a 
black-body of temperature [T^ 

r- W -(^)- 1 K-( i ^)- 1 G.v. (18, 

This assumes that the hole has no charge or an- 
gular momentum. This is a reasonable assumption 
since charge and angular momentum will also be lost 
through quant um e mission but on a shorter timescale 
than the mass |l2l| . This means that it loses mass at 
a rate 

M = -5 x 10 25 (M/g)- 2 f(M) g s" 1 (19) 

where the factor f(M) depends on the number of par- 
ticle species which are light enough to be emitted by 
a hole of mass M, so the lifetime is 

t(M) = 6 x 1(T 27 / (My 1 {M/gf s. (20) 

The factor / is normalized to be 1 for holes larger than 
10 17 g and such holes are only able to emit "mass- 
less" particles like photons, neutrinos and gravitons. 
Holes in the mass range 10 15 g < M < 10 17 g are 
also able to emit electrons, while those in the range 
10 14 g < M < 10 15 g emit muons which subsequently 
decay into electrons and neutrinos. The latter range 
includes, in particular, the critical mass for which r 
equals the age of the Universe. 

Once M falls below 10 14 g, a black hole can also 
begin to emit hadrons. However, hadrons are com- 
posite particles made up of quarks held together by 
gluons. For temperatures exceeding the QCD confine- 
ment scale of Aqcd = 250—300 GeV, one would there- 
fore expect these fundamental particles to be emitted 



rather than composite particles. Only pions would be 
light enough to be emitted below Aqcd- Since there 
are 12 quark degrees of freedom per flavour and 16 
gluon degrees of freedom, one would also expect the 
emission rate (i.e. the value of /) to increase dramat- 
ically once the QCD temperature is reached. 

The physics of quark and gluon emission from black 
holes is simplified by a number of factors [l04j . Firstly, 
one can show that the separation between succes- 
sively emitted particles is about 20 times their wave- 
length, which means that short range interactions be- 
tween them can be neglected. Secondly, the condition 
T > Aqcd implies that their separation is much less 
than Aq CD w 10~ 13 cm (the characteristic strong in- 
teraction range) and this means that the particles are 
also unaffected by strong interactions. The implica- 
tion of these three conditions is that one can regard 
the black hole as emitting quark and gluon jets of the 
kind produced in collider events. The jets will decay 
into hadrons over a distance which is always much 
larger than the size of the hole, so gravitational ef- 
fects can be neglected. The hadrons may then decay 
into astrophysically stable particles through weak and 
electomagnetic decays. 

To find the final spectra of stable particles emitted 
from a black hole, one must convolve the Hawking 
emission spectrum with the jet fragmentation func- 
tion. This gives the instantaneous emission spe ctrum 
shown in Fig.4 for a T = 1 GeV black hole [PO^ , The 
direct emission just corresponds to the small bumps 
on the right. All the particle spectra show a peak 
at 100 MeV due to pion decays; the electrons and 
neutrinos also have peaks at 1 MeV due to neutron 
decays. In order to determine the present day back- 
ground spectrum of particles generated by PBH evap- 
orations, one must first integrate over the lifetime of 
each hole of mass M and then over the PBH mass 
spectrum | l06j . In doing so, one must allow for the 
fact that smaller holes will evaporate at an earlier cos- 
mological epoch, so the particles they generate will be 
redshifted in energy by the present epoch. 

If the holes are uniformly distributed throughout 
the Universe, the background spectra should have the 
form indicated in Fig|5] All the spectra have rather 
similar shapes: an E~ 3 fall-off for E > 100 MeV due 
to the final phases of evaporation at the present epoch 
and an E- 1 tail for E < 100 MeV due to the frag- 
mentation of jets produced at the present and earlier 
epochs. Note that the E^ 1 tail generally masks any 
effect associated with the mass spectrum of smaller 
PBHs which evaporated at earlier epochs |27J. 

The situation is more complicated if the PBHs evap- 
orating at the present epoch are clustered inside our 
own Galactic halo (as is most likely). In this case, any 
charged particles emitted after the epoch of galaxy for- 
mation (i.e. from PBHs only somewhat smaller than 
M*) will have their flux enhanced relative to the pho- 
ton spectra by a factor £ which depends upon the halo 
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Log B (Kinetic Energy/GeV) 



FIG. 4: Instantaneous emission from a 1 GeV black hole 



concentration factor and the time for which particles 
are trapped inside the halo by the Galactic magnetic 
field. This time is rather uncertain and also energy- 
dependent. At 100 MeV one has £ ~ 10 3 for electrons 
or positrons and £ ~ 10 4 for protons and antiprotons. 
MacGibbon & Carr |l05j ] first used the observed cos- 
mic ray spectra to constrain f2pBH but their estimates 
have recently been updated. 
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FIG. 5: Spectrum of particles from uniformly distributed 
PBHs 



A. Gamma-Rays 

EGRET observations |l33fl give a 7-ray background 

of 

^1 = 7.3 x HT 14 ( 1 cm^GeV- 1 

dE \ lOOMeV J 

(21) 



between 30 MeV and 120 GeV. This leads to an upper 
limit H3 

^pbh < (5.1 ± 1.3) x 10~ 9 /i~ 2 , (22) 

where h is the Hubble parameter in units of 100. This 
is a refinement of the original Page-Hawking limit, but 
the form of the spectrum suggests that PBHs do not 
provide the dominant contribution. If PBHs are clus- 
tered inside our own Galactic halo, then there should 
also be a Galactic 7-ray background and, since this 
would be anisotropic, it should be separable from the 
extragalactic background. The ratio of the anisotropic 
to isotropic intensity depends on the Galactic longti- 
tude and latitude, the Galactic core radius and the 
halo flattening. Wright clai ms t hat such a halo back- 
ground has been detected [138J . His detailed fit to 
the EGRET data, subtracting various other known 
components, requires the PBH clustering factor to be 
(2 — 12) x 10 5 /i _1 , comparable to that expected. 



B. Antiprotons and Antideuterons 

Since the ratio of antiprotons to protons in cos- 
mic rays is less than 10~ 4 over the energy range 
100 MeV - 10 GeV, whereas PBHs should produce 
them in equal numbers, PBHs coul d only contribute 
appreciably to the antiprotons |135| . It is usually as- 
sumed that the observed antiproton cosmic rays are 
secondary particles, produced by spallation of the in- 
terstellar medium by primary cosmic rays. However, 
the spectrum of secondary antiprotons should show a 
steep cut-off at kinetic energies below 2 GeV, whereas 
the spectrum of PBH antiprotons should increase with 
decreasing energy down to 0.2 GeV. Also the antipro- 
ton fraction should tend to 0.5 at low energies, so these 
features provide a distinctive signature |95|. 

MacGibbon & Carr originally calculated the PBH 
density required to explain the interstellar antipro- 
ton flux at 1 GeV an d fou nd a value somewhat larger 
than the 7-ray limit [1.05] . More recent data on the 
antiproton flux belo w 0.5 GeV comes from the BESS 
balloon experiment |l42j and Maki et al. |lll| have 
tried to fit this data in the PBH scenario. They model 
the Galaxy as a cylindrical diffusing halo of diame- 
ter 40 kpc and thickness 4-8 kpc and then use Monte 
Carlo simulations of cosmic ray propagation. A com- 
parison with the data shows no positive evidence for 
PBHs (i.e. there is no tendency for the antiproton 
fraction to tend to 0.5 at low energies) but they re- 
quire the fraction of the local halo density in PBHs 
to be less than 3 x 10~ 8 and this is stronger than the 
7-ray background limit. A more recent attempt to fit 
the observed antiproton spectrum with PBH emission 
comes from Barrau et al. and is shown in FigEI 
PBHs might also be detected by their antideuteron 
flux 0. 
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PBH antiproton spectra 




antiproton kinetic energy (GeV) 



Some people have emphasized the possibility of de- 
tecting very high energ y co smic ray s from PBHs using 
air shower techniques |42l loo , l98l. However, this is 
refuted by the claim of Heckler [76| that QED inter- 
actions could produce an optically thick photosphere 
once the black hole temperature exceeds T cr .;t = 45 
GeV. In this case, the mean photon energy is reduced 
to m e (TBH /Tern) 1 / 2 , which is well below Tbh, so the 
number of high energy photons is much reduced. He 
has proposed that a similar effect may operate at even 
lower temperatures due to QCD effects 77]. Several 
groups have examined the implications of this pro- 
posal for PBH emission mil. However, these ar- 
guments should not be regarded as definitive since 
MacGibbon et al. claim th at Q ED and QCD inter- 
actions are never important |l08j . 



VI. PBHS AS A PROBE OF QUANTUM 
GRAVITY 



FIG. 6: Comparison of PBH emission and antiproton data 
from Barrau et al. 



C. PBH Explosions 

One of the most striking observational consequences 
of PBH evaporations would be their final explosive 
phase. However, in the standard particle physics pic- 
ture, where the number of elementary particle species 
never exceeds around 100, the likelihood of detecting 
such explosions is very low. Indeed, in this case, ob- 
servations only place an upper limit on the explosion 
rate of 5 x I0 8 pc" 3 y" 1 [1 Il29| . This compares to 
Wright's 7-ray halo limit of 0.3 pc~ 3 y~ 1 |13S| and the 
Maki et al. antiproton limit of 0.02 pc~ 3 y _1 

However, the physics at the QCD phase transition 
is still uncertain and the prospects of detecting explo- 
sions would be improved in less conventional particle 
physics models. For example, in a Hagedorn-type pic- 
ture, where the number of particle species exponen- 
tiates at the the quark-hadron temperature, the up- 
per limit is reduced to 0.05 pc _3 y _1 jSJ). Cline and 
colleagues have argued that one might expect the for- 
mation of a QCD fireball at this temperature |3lll3l| 
and this might even explain some of the short period 
7-ray bursts observed by B ATSE |4jj . They claim to 
have found 42 candidates of this kind and the fact 
that their distribution matches the spiral arms sug- 
gests that they are Galactic [3{J. Although this pro- 
posal is speculative and has been disputed [5^ , it has 
the attraction of making testable predictions (eg. the 
hardness ratio should increase as the duration of the 
burst decreases). A rather different way of producing 
a 7-ray burst is to assume that the outgoing charged 
particles form a plasma due to turbulent magnetic 
field effects at sufficiently high temperatures [lq . 



In "brane" versions of Kaluza-Klein theory, some 
of the extra dimensions can be much larger than 
the usual Planck length and this means that quan- 
tum gravity effects may become important at a much 
smaller energy scale than usual. If the internal space 
has n dimensions and a compact volume V n , then 
Newton's constant Gat is related to the higher di- 
mensional gravitational constant Go and the value 
of the modified Planck mass Mp is related to the 
usual 4-dimensional Planck mass M4 by the order- 
of-magnitude equations: 

G N ~ G D /V n , M™+ 2 ~ Ml/V n . (23) 

The same relationship applies if one has an infinite 
extra dimension but with a "warped" geometry, pro- 
vided one interprets V n as the "warped volume". In 
the standard model, V n ~ 1/M™ and so Mp ~ M4. 
However, with large extra dimensions, one has V n 3> 
l/M% and soM P «M 4 . 

A. Black Hole Production at Accelerators 

One exciting implication of these scenarios is that 
quantum gravitational effects may arise at the exper- 
imentally observable TeV scale. If so, this would have 
profound implications for black hole formation and 
evaporation since black holes could be generated in 
accelerator experiments, such as the Large Hadron 
Collider (LHC) . Two partons with centre-of-mass en- 
ergy y/s will form a black hole if they come within 
a distance corresponding to the Schwarzschild radius 
rs for a black hole whose mass Mpn is equivalent to 
that energy [ZfiL IHfj Il28j ] . Thus the cross-section for 
black hole production is 

<Jbh « 7rr|9(^ - M%%) (24) 
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where Mg^ is the mass below which the semi- 
classical approximation fails. Here the Schwarzschild 
radius itself depends upon the number of internal di- 
mensions: 



rs 



1 



M 



BH 



l/(l+n) 



(25) 



so that ubh °t s 2 /( n+1 \ This means that the cross- 
section for black hole production in scattering exper- 
iments goes well above the cross-section for the stan- 
dard model above a certain energy scale and in a way 
which depends on the number of extra dimensions. 

The evaporation of the black holes produced in this 
way will produce a characteristic signature |45l lo(l 
128] because the temperature and lifetime of the black 
holes depend on the number of internal dimensions: 



Tbh 



1 



rs 



tbh 



1 ( Mbh 
M P V M P 



(n+3)/(n+l) 



(26) 

Thus the temperature is decreased relative to the stan- 
dard 4-dimensional case and the lifetime is increased. 
The important qualitative effect is that a large frac- 
tion of the beam energy is converted into transverse 
energy, leading to large-multiplicity events with many 
more hard jets and leptons than would otherwise be 
expected. In principle, the formation and evaporation 
of black holes might be observed by the LHC when 
it turns on in 2007 and this might also allow one to 
experimentally probe the number of extra dimensions. 
On the other hand, this would also mean that scat- 
tering processes above the Planck scale could not be 
probed directly because they would be hidden behind 
a black hole event horizon. 

Similar effects could be evident in the interaction 
between high energy cosmic rays and atmospheric nu- 
cleons. Nearly horizontal cosmic ray neutrinos would 
lead to the production of black holes, whose decays 
could generate deeply penetrating showers with an 
electromagnetic component substantially larger than 
that expected with conventional neutrino interactions. 
Several authors have studied this in the context of the 
Pierre Auger experiment, with event r ates in excess 
of one per year being predicted [|| 1 128]. Indeed 
there is a small window of opportunity in which Auger 
might detect such events before the LHC. 



B. PBHs and Brane Cosmology 

The black holes produced in accelerators should not 
themselves be described as "primordial" since they do 
not form in the early Universe. On the other hand, 



it is clear that the theories which predict such pro- 
cesses will also have profound implications for PBHs 
because, at sufficiently early times, the effects of the 
extra dimensions must be cosmologically important. 
Such effects are not fully understood but studies of 
PBH formation and evaporation in brane cosmology 
are beginning to explore them. 

Brane cosmology modifies the standard PBH for- 
mation scenario. In particular, in Randall- Sundrum 
type II scenarios, a p 2 term appears in the Friedmann 
equation and this dominates the usual p term for 



T > 10 18 (Z/Z P )- 1/4 GeV 



(27) 



where I is the scale of the extra dimension. This 
exceeds 1 TeV for I < 0.2 mm. Black holes which 
are smaller than I are effectively 5-dimensional. This 
means that they are cooler and evaporate more slowly 
than in the standard scenario: 



Tbh oc M~ 1/2 , t cx M 2 



(28) 



The critical mass for PBHs evaporating at the present 
epoch is in the range 3 x 10 9 g to 4 x 10 14 g. This 
modifies all the standard PBH evaporation constraints 
[3(||65|,|8!| and, if some cosmic rays come from PBHs, 
it me ans that these could probe the extra dimensions 
130]. Another complication is that accretion could 
dominate evaporation during the era when the p 2 term 
dominates the p term in the Friedmann equation |l!3| . 



VII. CONCLUSIONS 

We have seen that PBHs could be used to study 
the early Universe, gravitational collapse, high en- 
ergy physics and quantum gravity. In the "early Uni- 
verse" context, useful constraints can be placed on 
inflationary scenarios and gravitational memory. In 
the "gravitational collapse" context, the existence of 
PBHs could provide a test of critical phenomena. In 
the "high energy physics" context, information may 
come from observing cosmic rays from evaporating 
PBHs since the constraints on the number of evaporat- 
ing PBHs imposed by gamma-ray background obser- 
vations do not exclude their making a significant con- 
tribution to the Galactic flux of electrons, positrons 
and antiprotons. Evaporating PBHs may also be de- 
tectable in their final explosive phase as gamma-ray 
bursts if suitable physics is invoked at the QCD phase 
transition. In the "quantum gravity" context, the for- 
mation and evaporation of small black holes could be 
observable in cosmic ray events and accelerator exper- 
iments if the quantum gravity scale is around a TeV. 
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